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The evolution of local spin polarization in finite systems 
involves interference phenomena that give rise to quantum 
dynamical echoes and non-ergodic behavior. We predict 
the conditions to observe these echoes by exploiting the NMR 
sequences devised by Zhang et al. [Phys. Rev. Lett. 69, 2149 
(1992)], which uses a rare 13 C as local probe for a dipolar 
coupled 1 H spin system. The non-ideality of this probe when 
testing mesoscopic systems is carefully analyzed revealing the 
origin of various striking experimental features. 

In recent years interest has aroused on transport prop- 
erties of mesoscopic systems [Q. Examples range from 
photons propagating through a colloidal solution to elec- 
trons in nanostructures. In most of the systems the quan- 
tum phenomena manifest as transport properties which 
are not linear with the size of the system. Since the 
Boltzmann equation does not retain interferences after 
collisions, it is necessary the use of the Schrodinger equa- 
tion (e.g. in Liouville or Keldysh's form Q). Thus, 
a localized excitation, after diffusing away and reaching 
the boundaries, might return as quantum dynamical 
echoes [Q resembling the familiar sound echoes. A triv- 
ial case occurs when a particle is placed in one side of a 
symmetrical double well potential and oscillates among 
them H . While in these examples the high group veloc- 
ities of the excitations restrict the observation of inter- 
ference phenomena in the time domain, the local excita- 
tions in a system of nuclear spins with magnetic dipolar 
interactions evolve with characteristic times on the hun- 
dreds of micro-seconds and can be tested by using an 
NMR technique devised recently by Zhang, Meier and 
Ernst Q (ZME). While they addressed a realization of 
the Lochsmidt daemon Q, they showed that, by applying 
an appropriate radio-frequency pulse sequence involving 
cross-polarization ||, the spin S of a rare 13 C (1.1% 
abundance) bonded to a with spin I± can be used as 
a local probe that injects magnetization in the proton 
and captures it at a later time. The set of magnetic dipo- 
lar coupled 1 H within a molecule constitutes the meso- 
scopic system where the spin dynamics can be monitored. 
The rest of the crystal constitutes a weakly interacting 
thermal bath. 

The experimental j^] time evolution of the polarization 
of a proton spin in a crystal of ferrocene, (Cstls^Fe, 
shows a nearly exponential decay of the local magneti- 
zation toward an asymptotic value of 0.2 of the initial 
polarization. This was interpreted as an "unhindered" 
spin diffusion among the five protons of a single ring. 
This is based on the assumption that the 13 C is an ideal 



local probe, i.e. the polarization transfer between the 
1 H and 13 C is complete and it takes a short time as com- 
pared with the characteristic time of the 1 H- 1 H interac- 
tions j|. However, there are two aspects in the results 
of this ingenious experiment which do not fit with our 
current understanding of mesoscopic quantum dynamics: 
One is the linear decay of the magnetization observed for 
short times which is obviously inconsistent with the ex- 
pected quadratic decay. The other is the monothonic de- 
cay which does not reveal any of the interference features 
characteristic of quantum evolution in systems with few 
degrees of freedom. The search for the conditions to ac- 
tually observe quantum dynamical echoe and answer to 
these questions lead us to solve the Schrodinger equation 
for the model system proposed by ZME. 

Let us consider an ideal situation in which there is 
a perfect local probe which is able to inject magnetiza- 
tion in a single -"-H nucleus and to measure it later. The 
system consists of N /-spins governed by the magnetic 
dipolar [a = 1) Hamiltonian truncated with respect to 
the dominant Zceman interaction: 
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Site j spin operators Ij are dimensionless. Subscripts 
start sequentially from J = 1. The interaction parameters 
are 



Mo 7? 



h 2 1 



- (3 cos 2 
2 x 



-1) 



(2) 



where 7 7 is the gyromagnetic factor of the /-spins, r^ 
are internuclear distances and 6jk their angles with the 
static magnetic field. A simplifying fact about ferrocene 
is that the ring performs fast rotations around its five- 
fold symmetry axis with a very short correlation time 
JlOj l (w 10~ 12 s) and therefore the interaction parameter 
is time averaged. The new constants depend fill only 
on the angle between the molecular axis and the mag- 
netic field and the angles jj t k between the internuclear 
vectors and the rotating axis. The factor in brackets in 
Eq.(|) becomes (|(3 cos 2 8 - 1)(3 cos 2 7j - k - 1)) = 1/2 
for a magnetic field normal to the axis. All sites in a 
molecule are then equivalent and the five fold symmetry 
is recovered, giving d\ — 1576Hz x 2irh for nearest neigh- 



bors. A locally polarized initial state 
a final state |/) with probability: 



(f | exp[--H//ti] 



evolves toward 



(3) 



2 



In the high temperature limit all the Ni = Nf configura- 
tions compatible with an up projection of the spin I\ are 
equally probable. Then the magnetization will be: 



M(t) = 2 
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For an ideal system of five spins arranged in a ring con- 
figuration, there is 1 configuration with total projection 
5/2, and there are 5, 10, 10, 5 and 1 configurations with 
projections 3/2, 1/2, —1/2, —3/2 and —5/2 respectively. 
On each of these subspaces there are 1, 4, 6, 4, 1 and 
possible initial and final states with spin I\ polarized 
summing up Ni — 16. The evolution under Hamiltonian 
(|l|) is called 'diffusion' in the laboratory frame. For the 
results of local polarization presented in Fig. 1 we multi- 
ply by — [\] the above interaction constants, which corre- 
spond to the experimental case of 'diffusion' in the rotat- 
ing frame. The dashed line represents the evolution when 
only nearest neighbor interactions are taken into account, 
while the full line considers all the interactions. At inter- 
mediate times, the first shows a clear quantum dynami- 
cal echo at around 1400^s « 5 x y/2 x h/([^]di), which 
can be interpreted as the time the excitation needs to 
wind around a ring of length L = 5a at an average speed 
«a//\/2, with maximum group velocity vm = a[^]di/h. 
This is confirmed by the study of correlation functions 
at different sites: the polarization at site 2 and 5 has 
the first maximum at 280/is and at sites 3 and 4 around 
560/is. Including second neighbors interactions, maxima 
are split and weakened, this can be assigned to the ef- 
fect of short-cutting the path around. Although dynam- 
ical echoes are still visible, one may conclude that they 
are weaker when more degrees of freedom become effec- 
tively mixed. In the long time regime a striking effect 
is noticeable in the inset of Fig. 1: local polarization 
does not distribute around the value 1/5 but rather os- 
cillates around M=lim t ^oo i" 1 f*M(t')dt' w 0.32. This 
non-ergodic behavior is a consequence of the degeneracy 
of opposite quasi-momenta in perfect rings. We have ob- 
served this for single rings as well as for pairs of coaxial 
rings with N < 10 spins interacting according to dipo- 
lar (a — 1), Ising like (a ^> 1), Heisenberg (a = — ~) 
and XY (a — 0) interactions by studying M. In the 
last situation we obtained the analytical high temper- 
ature solution of Eq. (0), which for odd rings gives: 
Mi = 1/N + (N — l)/N*, where the second term ac- 
counts for the degenerate quasi-momentum states. The 
complete ferrocene molecule gives M\ ~ 0.2. Here we dis- 
tinguish rings rotating independenly or as a whole, keep- 
ing an eclipsed or staggered conformation. While the first 
present a weaker version of the dynamical echoes of the 
single ring, in the other two the echoes and valleys merge 
in a sort of plateau for intermediate times. For very short 
times the magnetization decreases quadratically with t 1 . 
The coefficient is the configuration average of the sum of 
squares of the coupling constants effective to the flip-flop 



process in each configuration. 

These results clearly do not represent the experimental 
data in ZME Although quantum dynamical echoes 
are indeed attenuated by including the interactions with 
other rings and molecules, the linear decay of the mag- 
netization suggests that some quantum evolution has al- 
ready occurred for the experimental time ti = 0. There- 
fore, it is important to consider that magnetization is 
injected and monitored using the spin S, which consti- 
tutes a non-ideal probe. The complete pulse sequence 
is schematized in Fig. 2. The most relevant part is as 
follows: A) The magnetization from an initially y po- 
larized S spin is transferred to the y-axis of the I Y spin 
(and in lesser degree to its neighbors) through a cross po- 
larization pulse of duration td, when both the abundant 
and rare spins are irradiated at their respective resonant 
frequencies with field strengths fulfilling the Hartmann- 
Hahn || condition JiBfj = jsB^ s - The time td is the 
shortest of those that maximize the polarization transfer 
M. B) The /-spins evolve in presence of B v u during a 
time t\. Thus, the relevant evolution of the spin-diffusion 
sequence occurs in a rotating frame. C) Another cross- 
polarization pulse of length t p = td is applied to transfer 
back the polarization to the a;-axis of S. D) The S polar- 
ization is detected while the /-system is kept irradiated 
(high-resolution condition). Therefore a better simula- 
tion of the experiment includes the /'s spins system, the 
non ideal probe S and their interaction. The IS dipolar 
interaction extends the sum in (§) to the S spin. The in- 
teraction coefficient is b k defined as (^) with a ji replaced 
by a 7s and distances and angles modified accordingly. 
Given the difference in gyromagnetic factors the dipo- 
lar IS interaction is truncated dropping terms that mix 
different Zeeman subspaces, hence: 
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When two r.f. fields are applied to both / and S 
spins, these are quantized in their corresponding rotating 
frames and the total Hamiltonian is: 



H y{x) =Hn+Hi+n 
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the superscript indicates the direction of B 1S . 

The isolated I-spins system is described by the first 
two terms, where Hi = —Tmju J2k rfl with hu>u = ~iiB v XI 
is the effect of the r.f.. In the case of strong irradiation, 
| fawn |3>| djk | , we can approximate (nl) in a frame 
tilted to the rotating y axis by 
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As usual, we neglected non-secular terms of the form 
It is the factor — [~] relating (Q) and 
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what allows the realization of the Lochsmidt dae- 
mon (?| by experimentally switching the evolution from 
one situation to the other which leads to a polarization 
echo B, not to be confused with the quantum dynami- 
cal echoes described above. The isolated probe in the 
presence of a r.f. field along y (x) direction is described 

by H y J x) = -hui S S y W with haj ls = jsB^g\ The 
system-probe interaction with an additional trunca- 
tion in the doubly rotating frame is: 

n IS « 5>f [i£s- + (-1) / fe -s+], (8) 

k 

where factors in parenthesis correspond to Bf s . In Fig. 
2, the evolution starting with the polarized S spin in- 
volves two Hamiltonians, during the times td and t p the 
evolution is given by the total Hamiltonian Tt yt - X \ but 
during t\ as well as during the acquisition, the irradia- 
tion makes the action of Tt IS negligible (high resolution 
condition) and evolves with r H! II + Hj . The probability 
of having an initial state \i v ) with the probe polarized 
along y direction and a final state \f x ) with the probe 
polarized along x is: 

PfM =1 (fx I cxp[-iW%]exphi(W, + H' II )t 1 ]exp[- 

(9) 

The resulting magnetization is shown in Fig. 3 with 
open circles. For short t\ the evolution is linear with time 
and the dynamical echoes are still noticeable but appear 
shifted to shorter times. Both effects are consequence 
of the evolution of the magnetization in the I-I system 
during the two cross polarization periods: by comparing 
with Fig.l, it is clear that the overall effect of them is 
to shift the curve in a time shorter than their sum. This 
is verified by superposing the ideal evolution from Fig. 
1, shifted in £ s hift = 126/is < (td + tp) and renormalized 
with a factor of a — 1.22 obtained from the short time 
behavior. This explains why the asymptotic value of the 
magnetization in the complete sequence has been raised 
to M n .i. = oMj m 0.4 from the Mi w 0.32 around which 
the ideal polarization of Fig.l oscillates. Then M n .i over- 
estimates the fraction of magnetization remaining in the 
proton. The fact that Figs. 1 and 3 can be superposed 
indicates that S-I interactions do not break the five-fold 
symmetry appreciably within the relevant experimental 
times, this is a consequence of the high resolution condi- 
tion. 

These results indicate that the experimental H asymp- 
totic value M cxp = 0.2 = M n .i originates on the evo- 
lution that has already occurred at t\ = 0,when po- 
larization starts to be recorded and is normalized. A 
more proper normalization would give an ideal asymp- 
totic value M\ w M cxp /a = 0.16, lower than M\ for the 
complete molecule (two coupled rings). Since, as dis- 
cussed above, the quantum evolution does not give an 
equi-distribution of polarization, to achieve the ex- 
perimental 0.2, the evolution must involve a system with 



TV 3> 10. With the participation of an important num- 
ber of spins other than those belonging to one molecule, 
the boundary of the system becomes fuzzy, weakening 
the interferences that originate the quantum dynamical 
echoes. Then, it becomes clear why these do not show 
up in the reported experiments. 

An interesting phenomenon shown up by the non-ideal 
evolution in Fig 3 is the small high frequency (w 
oscillation around the superposed ideal curve. This ef- 
fect is due to an incomplete S-I transfer which allows 
interference of phases accumulated during evolution t\ 
by sub-spaces with different total spin ^2 I k y projection; 
i.e. the interference of the dashed and dotted paths of po- 
larization amplitudes in Fig. 2. They may originate the 
small 'dispersion' observed in the top of the polarization 
echo in the experimental curve. 

In summary, we have solved the relevant part of a non- 
ideal spin 'diffusion' experiment for a simple five-spins 
system resorting to the usual assumptions to generate 
the system, probe and system-probe Hamiltonians: trun- 
cation of the dipolar Hamiltonian and total decoupling 
under high resolution conditions (during the time t\). We 
have shown that while 13 C behaves as a non-ideal local 
jarobe, its 2 effect can be reasonably quantified. As occurs 
S^wiiJllI naesDscopic physics the probes are actually part of 
the system |L2] and therefore have a relevant role in the 
interpretation of the results. This undesired effect can be 
experimentally overcome in part, by studying the diffu- 
sion in the laboratory frame which tends to compensate 
the non- ideality of the probe and shrinks in a factor [i] 
the time scale, 'cleaning' the evolution from interactions 
different than the secular dipolar one. We predict that 
a spin 'diffusion' experiment in molecular crystals con- 
taining few interacting spins per molecule should present 
quantum dynamical echoes. Possible candidates are 
molecules with single rings of the form (C„H„). Molecules 
with linear topology in the interaction network present 
less important second neighbors interaction, making the 
interference effects even stronger. 

This work was performed at LANAIS de RMN (UNC- 
CONICET) with financial support from Fundacion An- 
torchas, CONICOR and SeCyT-UNC. 
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Figure 1: Ideal evolution of a 1 H-spin polarization 
function of t\ for a system of 5 spins in a ring. 
Dashed line considers only nearest neighbor interaction 
d\. Full line is the exact dipolar evolution. The inset 
shows fluctuations around the value M\ = 0.32. 



Figure 2: Pulse sequence for measurement of proton 
spin 'diffusion' in the rotating frame. An/2 pulse on the 
abundant X H spins system creates a polarization that is 
transferred during tc to a rare 13 C. After the decay of 
the proton spin coherence during ts- (A) the 13 C po- 
larization is locally injected to the bonded proton, (B) 
'diffusion' in the protons system is allowed and (C) it is 
captured again in the 13 C where (D) it is recorded. 



Figure 3: Non-ideal evolution of a 1 H spin polariza- 
tion in a ring of ferrocene (open circles) as detected in 
a strongly coupled 13 C spin. Calculation uses the pulse 
sequence in Fig. 2 with td — t p — 85/zs with the molec- 
ular symmetry axis at an angle of n/2 with the external 
magnetic field. Full line is the ideal evolution in Fig.l, 
shifted in 126/xs and normalized accordingly. 
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